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We analyze the optimal state under the canonical phase measurement, as given by Berry and 
Wiseman [Phys. Rev. Lett 85, 5098, (2000)], in the presence of photon loss. The model of photon 
loss is a generic fictitious beam splitter, and we present the full density matrix calculations, which 
are more direct and do not involve any approximations. We find for a given amount of loss the 
upper bound for the input photon number that yields a sub-shot noise estimate. 
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I. INTRODUCTION 

Canonical phase measurement in quantum mechanics 
is a significant problem, for the main reason that phase is 
a quantity that is conjugate to the number, TV, of photons 
in a particular electromagnetic mode [TJ HJ |3] . Due to this 
conjugate nature, the phase estimate Aip is ultimately 
limited by the number N of the photons as A(p — 1/N, 
which is conventionally referred as the Heisenberg limit. 
In the usual classical setting, such as interferometry with 
lasers, for a given number of input resources, TV, phase 
estimate scales as 1/^/N, which is usually referred as 
shot-noise limit. 

Accurate phase estimation has many practical applica- 
tions such as metrology, imaging and sensing 14 . Achiev- 
ing Heisenberg limit in practice is not a trivial problem 
and there have been numerous proposals to achieve this 
limit [51[i|71iaini[TnilIIllIlllIl[Tl[Tn]. canonical phase 
measurement has been first dealt with Helstrom [TB] and 
Shapiro [17 , and later Sanders and Milburn [HI [19] used 
it to obtain a phase estimate in a Mach-Zehender inter- 
ferometer (MZI) as shown in Figjl] (excluding the loss 
part). The phase estimate thus obtained is independent 
of the system phase, unlike in other methods [SI [3 |7] 
where the ultimate limit is achieved for a particular sys- 
tem phase. Also, the measurement specified by Sanders 
and Milburn is not particular to a specific input state. 

Motivated by this work. Berry and Wiseman [8 an- 
alytically derived an input state, called as the optimal 
state, subject to the canonical measurement. They also 
suggested an adaptive method of approximately imple- 
menting the canonical phase measurement [SI [50] , which 
has been used in the recent experimental realization of 
the Heisenberg limited phase measurement by Higgins et 
al. |21j . The canonical measurement, written as a Pos- 
itive Operator Valued Measure (POVM) by Sanders and 
Milburn [TH] is, 

2j + 1 

F{ip)d(p = — IjV) (jVI dip, (1) 

in terms of the phase states 

a=-j 



In defining this, the Schwinger's representation is used, 
and for completeness we wish to outline the notation. 
The three angular momentum components Jx , Jy and Jz 
are very effective in analyzing two-port, lossless, interfer- 
ometers [5J[B]|TS]. For the two modes, d and b of the MZI 
(Fig. [ij, these two mode operators are, 

X = {d^b+dP)/2, Jy = {a^b-aP)/2i, 
J, = (fltfl - Pb)/2, J^^J^ + J^ + Jl 

In the context of the MZI, implements the opera- 
tion of a 50-50 beam splitter as e^i'^/^)^'' and Jz defines 
the photon number difference in two modes. The simul- 
taneous eigenvector of and Jz, \j,'m)z represents the 
joint input state \j + m)^, and \j — rn)i^, in the Fock- 
state basis, and the total input number of photons is 
N = 2j. The simultaneous eigenvector of and Jy, 
which is \j,n)y, represents the joint state within the in- 
terferometer. The beamsplitter transformation in this 
representation performs a rotation about the Jx- The 
phase states discussed above, are defined in terms of 
states within the interferometer in Fig. [Tj and thus the 
output modes or detectors are irrelevant for the present 
purpose. Thus the probability distribution for the system 
phase (j) is obtained as 

P{^)d^ = (V'l F{ip) IV) dip ^ Tr[pF{^)]dip. (2) 

Note that ip is the estimate of the system phase (p. The 
optimal state, to be specified below, is derived condi- 
tioned upon minimizing the Holevo variance calculated 
from the above probability distribution. The Holevo 
phase variance is defined as [22 

{Ap)'^-l + \(e^n\-^ (3) 

where |(e*^)| = J^'' dipP{ip)e''^-f , is also called the 
sharpness. Here ip is a, mean phase and we take it to 
be zero. The Holevo variance for the optimal state is 
given by Ref. fElW\, 

(Av»^^ta„^(^)»^, (4) 

thus giving rise to a phase estimate that scales as the 
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In this paper we study the performance of optimal 
state and the optimal POVM in the presence of the pho- 
ton loss associated with propagation. We use the generic 
beam splitter model for photon loss [23", as shown in 
Fig. [l] The input mode c for this fictitious beam splitter 
is a vacuum mode, and the output mode c' is then to 
be traced out. This typically implies that the photons 
that are lost in mode c', due to the nonzero reflection 
coefficient r, correspond to the photon loss. 

In the following two Sections we describe optimal state 
in presence of photon loss and carry out the explicit den- 
sity matrix calculation. In Section [IV| we quantitatively 
describe the effect of photon loss on the canonical phase 
estimation. Section |V] concludes with numerical results 
and discussions. 



Loss 



FIG. 1: Schematic of Mach-Zehender Interferometer (MZI) 
with photon loss at the phase shift. The state \tp)in represents 
the joint input state at ain and bin- Photons in the lower arm 
first encounter the fictitious beam splitter (BS), for which 
vacuum enters through the other input, and depending the 
transmission coeflRcient, some of them are scattered into the 
mode c', which are ignored (traced out), and the remaining 
pass through the phase shifter. 



Heisenberg limit. Besides the optimal state, other promi- 
nent states that achieve the Heisenberg limit are so called 
the NOON states 0113] and the dual-fock state [71110]. 
Formally the NOON state is, 



/N 



\m,\oh + \o),\Nk 

V2 



(5) 



Note that the above state is not an input state at the MZI 
shown in Fig.l, but a state within the interferometer and 
the subscripts a and b denotes the internal modes of the 
MZI. The dual-fock input state is given as. 



(6) 



The detection scheme for both, the NOON and the 
dual-fock state that results in the Heisenberg limited 
phase estimate is the parity detection pUl , first pro- 
posed by Bollinger et al. [T3] in the context of frequency 
metrology with trapped ions. The shot- noise and sub 
shot-noise limit with matter wave interferometry which 
deals with Bose-Einstein condensates at the input of the 
MZI is studied in Refs. [271. 

It is natural to question the performance of such states 
or the detection schemes in a more realistic conditions 
such as photon loss associated with the propagation. The 
analysis of the NOON states under propagation loss was 
carried out independently by Gilbert et al. |25j and by 
Rubin and Kaushik [26], where they used pure state for- 
malism. In Ref. [25] |26] for NOON state, the minimum 
number of photons required to achieve a minimum de- 
tectable phase in presence of loss is also given. 



II. OPTIMAL STATE IN PRESENCE OF 
PHOTON LOSS 



We now proceed to develop the mathematical frame- 
work to study optimal-state canonical interferometry 
with photon loss. The beam splitter representing loss, 
with arbitrary transmission and reflection, can be char- 
acterized by an angle 9, such that transmission and reflec- 
tion coefficients are r — cos^(6'/2) and r = 1 — cos^(6'/2), 
respectively. Therefore, the loss is simply the reflection 
coefficient. 



L = r= 1 -cos2(6'/2). 



(7) 



The action of such an arbitrary beam splitter on an ar- 
bitrary joint input state | j, a) in Schwinger notation, is 
simply given as [B], 

e'«^Mj,a)= ^ e^''^'-^^dl,{e)\j,b), (8) 

b=-j 

where ^(^) is the usual rotational matrix element given 
as: 



dUiO) = (-1) 



a—br 



x(l - COS 6*) 2 (1 + cos 



where Pi"'^\x) is the Jacobi polynomial [28j. Also 
it is worth noting that when converted to Fock-basis, 
a two mode joint state in Schwinger notation, |j, a) is 
\j + a) \ j-a). 

The optimal state was originally derived by Berry and 
Wiseman [8] conditioned on minimizing the phase vari- 
ance with the canonical probability distribution given in 
Eq. ([2|. Formally the optimal state is. 



I {j - ay.jj + ay. (a-b,a+h) , 



(9) 



opt 



E 



2j + 2 



b^).- (10) 
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Recall that the simultaneous eigenstate of and jy de- 
note the state within the interferometer, and thus the 
above state is after the first beam splitter, while \j + IJ) ^ 
and \j — fi)f^ represents the fock state corresponding to 
modes a and b respectively. Rewriting the above state in 
a more explicit form as the product of the states at the 
two arms of the interferometer as, 



/ opt 



where 



b ' 



(11) 



2j 



With the loss in mode a, which is represented by the fic- 
titious beam splitter, and |0)^ is the state entering the 
other input port, the combined input state for the ficti- 
tious beam splitter is: |j + |0)^, and thus the state to 
be considered is. 



(12) 



The fictitious beam-splitter transforms modes a and c 
to modes a' and c' respectively. Making use of the 
Schwinger representation for modes a and c, the input 
b + f^a l^)c fictitious beam splitter can be writ- 

ten as a joint state: |^^, ^-^)ac Letting k = {j + fx)/2 
and using Eq. (|8| we have the output as, 

m— — /c 

k 

^ J2 e^i^"^-'Ui,{9)\k + m)^,\k-m)^,. (13) 

m— — k 

Therefore the pure state of the inner modes a', 5 of the 
interferometer, and mode c' of the lost photons is given 
as, 

m— — k 

X |fc + m)„, |fc-m)^, (14) 

where fc = (j + /i)/2 and d^^{6) is the usual rotational 
matrix element, as defined earlier. 



III. 



DENSITY MATRIX DESCRIPTION 



The state specified in Eq. (14 1 is a pure state and cannot be used further, so we need to calculate the reduced 



density matrix, by tracing out, mode c', as we have no more access to the lost photons. Thus first we need to calculate 
the total density matrix, representing the pure state of Eq. ( 14 1 , 



= E E E i^,^udl^mdiAO)e''^^"'-'''^e-^^^^--'^^^^ (15) 

f-i^u——j m— — k n— — k' 



where k' = {j + J^)/2. The total density matrix given in Eq. (15) explicitly represents the state within the 
interferometer for a given loss, characterized by the angle 9, and is useful in analyzing lossy interferometers such 
as the present one. 

As the mode c' is to be ignored, we need the reduced density matrix by tracing out that mode from the total density 
matrix given in Eq. (15). Thus we have, 

p' = Tr,4|^) (V'l] 

j k k' 

= E E E ^M^X,fc(e)<fc'(^))e^^^""''e-^^('-'=')|ft + m)^,b--Ai),,,(fc'-n|,(j-H^ (16) 



Noting that c' {k — m\k' — n)^' — 5k-m,k' ~m which eliminates the two exponential terms in Eq. ( 16 1. This leads to: 

j k k' 

P' ^ E E V'MV'.<fc(^)<fc'(^)4-™.fc'-n|fc + m)^, 

^M——j rn — ~k n——k' 



We can now use Eq. (17) in Eq. (|2| to obtain the probability distribution and thus the minimum detectable phase as 
a function of 0-which characterizes the photon loss-and the input photon number 2j. 



(17) 
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FIG. 2: Log-Log graph of minimum detectable phase versus the input photon number A'^ = 2j for two different values of loss, 
(a) L = 10~^ and (b) L = 0.3. The dotted line in both figures is the shot noise limit 1/^/2] and the dashed line in (a) is 
the Ifeisenberg limit given be Eq. Q, for comparison. The solid lines are numerically evaluated using Eq. ([3| and Eq. (24 1. 
Depending on the loss, the minimum detectable phase starts to diverge at certain photon numbers, which we shall call the 
optimal number, Nopt- 



IV. PHASE ESTIMATE IN PRESENCE OF 
PHOTON LOSS 

To get an estimate using Holevo phase variance, we 
need to calculate the probability distribution, P{ip) using 
Eq.([2|. We thus, in presence of photon loss have, 

P{^)=Tr[p'F{^)]. (18) 

Note that (p is an estimate of the system phase (p. 

The POVM, F{ip)dip [Eq. ^] is given in terms of the 
Jy eigen states. Noting that the Jy eigen states are the 
states of the modes within the interferometer, we can 
rewrite the POVM as. 



1 

2tt 



J2 e'^"-''^nj,a),(j,/?l 



Since the joint state \ j, a)y is written as \ j + a)^, \ j 
we have. 



(19) 



F{^) 



1 

2tj: 



3^-/3)^ \j-a). 



E 

®a' {j + PI {j - P\ (20) 

After the fictitious beamsplitter which represents pho- 
ton loss associated with propagation, the inner modes of 
the MZI are a' and b, and so the Jy eigen states are the 
product states of these modes. Thus using Eq. (20) in 
Eq. (18 1 along with Eq. (17 1 for the reduced density ma- 
trix, after carrying out the trace operation for the modes 
a' and b, we have the probability distribution in presence 
of loss as, 

j k k' 

^(^) = i E E E 



/i,^— — J 7n— — k n— — k' 



xe*("-^)'^(fc' - n\j + a){j + /3|fc + m) 

x{j-(3\j-fi){j-,y\j-a)]{2l) 



Recalling k' — {j + h')/2, k — {j + fJ,)/2, and because 
/3 = ^ and a = v as a consequence of the last two inner 
products in the above equation, we get from the first two 
inner products n = k' and m = k respectively. Hence the 
probability distribution explicitly given as. 



^(^) = ^ E ^^.V^.e^'^-^^^dtW^.^'W- (22) 



With this probability distribution we can calculate any 
moments of the phase estimate as a function of loss, 9. 
Here we calculate the Holevo phase variance rather than 
the standard variance, because the optimal state is de- 
rived by minimizing the Holevo phase variance [5] . Using 
Eq. (22) we now calculate sharpness Ke*"^)] is given as, 



2tc 



(23) 



where (5^^^_i = ^ " e(''-^+i)^d(/7. Also using 
the rotational matrix element in Eq. (|9]), we have 

4,ki^)dk',k'i^) = [cos2 f ] ^'"^'''^ and invoking the defi- 
nition of k and fc' in Eq. (23), we obtain 



(24) 



for the expression of the so-called sharpness as a function 
of loss. We use Eq. (24) in Eq. ([s]) to obtain the phase 
estimate as a function of loss. 
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FIG. 3: Log-Log plot of the optimal number, A'opt as a func- 
tion of loss for various values of loss L. The sub-shot noise 
range indicates the range of loss, so as to obtain the sub-shot 
noise phase estimate 



V. NUMERICAL RESULTS AND DISCUSSION 

Now the minimum detectable phase shift can be found 
with the uncertainty in phase estimate by plugging 
Eq. p4|) in Eq. In Fig. |2] we numerically plot the 



minimum detectable phase Atp as a function of 2j = N, 
the number of input photons, for two different values of 
loss, L = 10~^ and L = 0.3. For comparison we also 
plot the Heisenberg limit for the optimal state given in 
Eq. ^. We can immediately observe that the minimum 
detectable phase is not always a continuously decreasing 
function, depending on the loss, the minimum detectable 
phase starts to diverge at certain photon numbers, which 
we shall call the optimal number, N^pt. 

As we do not have a closed form expression for the 
Holevo variance, as a function of loss L and input photon 
number A'^, we have not found an analytical form for the 
optimal number, A^opt as function of loss, L [see Eq.(|7])] . 
So we numerically plot Aopt as a function of L in Fig. [3] 



and it clearly decreases as the loss L increases. 

Further, our analysis shows for small loss, about 0.05% 
loss, which is equivalent to 99.5% of transmission, we can 
still have the input number of photons, which gives a sub- 
shot noise limited phase estimate. 

On the other hand, the upper bound for a certain num- 
ber of input photons can be found where Aip meets the 
shot-noise limit in Fig. [2|a). That is to say, for a small 
amount of loss, we can still operate the interferometer, 
with the input number of photons less than the upper 
bound and achieve a sub-shot noise limited phase esti- 
mate. 

To summarize, we analyzed the canonical phase mea- 
surement in presence of photon loss. Our formalism is 
based on the density matrix, which describes the mixed 
states, which naturally arise due to the presence of loss. 
Our analysis shows that the minimum detectable phase is 
not monotonically decreasing function and would tend to 
increase at certain photon numbers, depending the loss 
present. Nevertheless, we can for small loss have con- 
siderably high number of input photons for the optimal 
state, that achieve sub-shot noise level phase estimates. 
We also have numerically plotted the optimal number, 
where the minimum detectable phase shift starts to di- 
verge, as a function of loss, L, which would determine 
the optimal number of the input photon number for the 
optimal state, in presence of given amount of loss. 
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